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Abstract 

We theoretically study the spin and charge generation along with the electron transport on a 
disordered surface of a doped three-dimensional topological insulator/magnetic insulator junction 
by using Green’s function techniques. We find that the spin and charge current are induced by not 
only local but also by nonlocal magnetization dynamics through nonmagnetic impurity scattering 
on the disordered surface of the doped topological insulator. We also clarify that the spin current 
as well as charge density are induced by spatially inhomogeneous magnetization dynamics, and 
the spin current diffusively propagates on the disordered surface. Using these results, we discuss 
both local and nonlocal spin torques before and after the spin and spin current generation on the 
surface, and provide a procedure to detect the spin current. 
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I. INTRODUCTION 


In spintronics, the mutual control of the direction and the flow of spin is a central issue 
for wide applications. The flow of spin, i.e., spin current, is the difference between the 
currents of up and down-spin conduction electrons. It is known that the spin current is 
induced in the setup of the ferromagnetic metal (FM)/normal metal (NM) junction- -3 . Its 
origin is due to the magnetization dynamics in ferromagnet, which transfers the spin angular 
momentum of the magnetization into that of the conduction electrons. The transfer of the 
spin angular momentum is called spin-pumping. Here, the spin-pumping of magnetization 
dynamics generates the spin-current in the NM, and the spin-current can be converted into 
charge current through spin-orbit interactions^ - -. 

Topological insulator (TI) is a new class of materials which has a gapless surface state, 
dubbed as the helical surface state, in which the spin and momentum are locked by the spin- 
orbit interactions^ - -. On the surface of the TI, the direction of charge current and that of 
the spin of conduction electrons can be mutually manipulated by an applied electromagnetic 
held through the spin-momentum locking. There have been many theoretical works in 
hybrid systems including superconducting junctions on the surface of TI stimulated by the 
exotic surface stated - — . In the TI/FM junctions, the anomalous charge-spin transport^ 2- — , 
the anomalous tunnel conductance^ - —, the giant magneto resistance--- - —, and the current- 
induced spin-transfer torque 24 ^ have been studied up to now. The exotic phenomena are 
triggered in the presence of static magnetization and an applied electromagnetic held. The 
magnetization of the ferromagnet plays the role of an effective vector potential for conduction 
electrons, which is like a vector potential of electromagnetic helds. Owing to the effective 
vector potential, the time-derivative of the magnetization can be regarded as an effective 
electric held, and the magnetization dynamics generates charge current on the surface of the 
TI/FM junction even in the absence of electromagnetic fields^ - — . This is called as the spin- 
charge conversion. The direction of the induced charge current is perfectly perpendicular 
to the magnetization dynamics due to the spin-momentum locking. The relation between 
the direction of the magnetization and the induced charge current can be a characteristic 
property on the surface of the TI. The property of the spin-pumping on the surface of TI 
can be applicable for spintronics devices. 

Existing works of the spin-charge conversion have been done in the case of a clean sur- 
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face of the TI, namely the ballistic transport regime. However, the actual charge trans¬ 
port on the surface of the TI is in the diffusive regime due to the nonmagnetic impurity 
scattering^ - —" . Since Burkov et al., have predicted not only the local but also the non¬ 
local current on the disordered surface of the TI in the presence of the applied electric field 22 ; 
we can naturally expect nonlocal current is driven by the magnetization dynamics even on 
the disordered surface of the TI in the presence of the magnetization dynamics. 

In this article, we study the charge-spin transport due to the magnetization dynamics on 
the disordered surface of the three-dimensional doped Tl/magnetic insulator (MI) junction, 
as shown in Fig. 1, where we show that charge current and spin polarization on the surface 
of the TI are induced not only by a local, but also by a nonlocal magnetization dynamics. 
Besides this, we clarify that the spin current is driven by the dynamics of the spatially 
inhomogenous magnetization, and the spin current diffusively propagates on the surface. 
The magnitude of the spin current reflects the spatially inhomogenous spin structure of the 
MI. The directions of the spin flow and the spin projection of the spin current are perfectly 
linked by the spin-momentum locking on the surface of TIs. The present features may serve 
as a guide to fabricate future spintronics devices based on the surface of TIs with magnetic 
substance. 

The merit of the choice of MI on the surface of the TI instead of metallic ferromagnet is 
to prevent the induced charge current going through the bulk of the MI. Then, we can focus 
on the charge transport on the surface of TI. Besides, the Gilbert damping constant in Mis 
tends to be smaller than that in ferromagnetic metals. The small value of the damping in 
Mis can be useful for the detection of the spin current on the disordered surface of the TI, 
as discussed in sec. IV Cl 

II. MODEL 

We consider conduction electrons coupled to an effective localized spin on the disordered 
surface of the three-dimensional doped TI attached with the MI, as shown in Fig. 1. The 
setup in Fig. 1 is similar to a system, where conduction electrons couple with the magnetic 
moments of ferromagnetic metals deposited on the surface of the TI—We expect 
that on the surface of the TI, the effective localized spin (S) can be produced from the 
magnetization in the MI through magnetic proximity effects. In the following, we use the 
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Hamiltonian, describing the surface of the TI with MI, given by 


U = 'HtI + U a d + Vimp, (1) 

where the first term in Eq. (|T|) , Uti is Hamiltonian of the conduction electrons on one of 
the surface of the doped TI without S as 

"Hti = j dx^{—ihv F (<r x V) z — ep]^, (2) 

Here, ^ = ^(x,t) = (ipL -0j), and ifj are the creation and annihilation operators of the 
conduction electron, respectively (where indices T and J, represent spin), ep is the Fermi 
energy, and is the Fermi velocity of the bare electron on the surface of the doped TI. The 
<r is the Pauli matrices in spin space. The second term of Eq. (IT]) ., U s d-> shows the exchange 
interaction between the conduction electron spin s = and the localized spin S on 

the disordered surface of the doped TI, as described by 

Usd = ~ J dxJshifts • crip, (3) 

where J s d > 0 is the exchange coupling constant. The localized spin S can be described 
by the magnetization of the MI as S — —(S/M)M, where S and M are the magnitude of 
the localized spin and of the magnetization, respectively. We consider that in general, the 
localized spin S = S(x,t ) depends on the time and position on the surface of the TI. The 
S(x,t) changes slowly compared with the electron transport relaxation time (r) and varies 
in space compared with the electron mean-free path (£). We expect that from the Eqs. 
o-o, the in-plane component of the localized spin, = S — S z z, can play the role of 
the effective vector potential for the conduction electrons on the surface. The out-of plane 
component of the localized spin S z plays a role to open the energy gap of the dispersion on 
the surface of the doped TI. We assume that the band gap opened by S z is smaller than the 
Fermi energy on the surface of the doped TI, he., ep — J s dS z > 0. The third term of Eq. (P), 

Ni 

3 = 1 

represents nonmagnetic impurity scattering on the disordered surface of the doped TI. The 
impurity scattering causes the relaxation time r of the transport of conduction electrons on 
the surface of the TI. Here U[ = U\8(x—r 3 ) is a delta-function type potential, u- x is a potential 


dxU^ip, 


( 4 ) 
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FIG. 1: (Color Online) A setup of spin-charge current generation due to magnetization dynamics 
in Mis deposited on disordered surface of a doped TI. 

energy density, r 3 is the position of impurities, and N\ shows the number of impurities. The 
contribution of V\ can be treated by the impurity average. 

We will calculate the spin current and charge current due to the spin-pumping in the 
linear response to S under the condition J sd <C h/r. We expect that the condition could 
be realized on a metallic disordered surface of the TI, which satisfies h/(e F r) -C 1. For 
example, the exchange coupling J sd can be estimated by J sd ~ 6meV— for Ni 81 Fei 9 /Bi 2 Te 3 
junction. If 10 2 is satished, the perturbation can be accessible. 

A. Renormalization of the Fermi velocity 

Green’s function on disordered surface of doped TI can be described by using TLt\ in¬ 
cluding Vjmp within the self-consistent Born approximation of bj mp as 


gk,uj = [huj - {hv F (& x k) z - e F } - E fcjU; ] \ 


(5) 


where T,k yU1 is the self-energy within the Born approximation given by 



( 6 ) 


k' 


Here, E*.^ satishes the Ward-Takahashi identity—. To estimate the value of Sfc jLJ , we consider 
the /^-dependence of Uk-k which plays the role to prevent the ultraviolet-divergence over 
a large momentum in the /^-integratio n 25 ' 34 . When can be described by 2 x 2 matrix 
= So + E" + E z a z , E is estimated, where E 0 and E z are independent of k, while 
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EH = E x a x + T, y a y depends on k. Then, the Green’s function ^ is described by E within 
the Born approximation as^— 


9k,u 


zH> 

hhj — x k) z — ep} H- 

2 T 


-i 


(7) 


where g r ku represents the retarded Green’s function. From Eq. (JTJ), the Fermi velocity is 
renormalized by E^, and the renormalized Fermi velocity is represented by vp = Vp/(1 + £), 
where £ = niU^/(4rch 2 Vp) is a small value depending on the relaxation time^. The last term 
in Eq. (El) is caused by the retarded component of Im[E 0 ], 

Equation (ED indicates the Green’s function on the disordered surface of the doped TI 
estimated within the Born approximation of V irav . Therefore, we could expect that hvp(& x 
k) z — ep in the Eq. (ED corresponds to the dispersion on the disordered surface of the TI. 
The dispersion is different from that of Hti without Idmp- In the following work, we will use 
an effective Hamiltonian "Hti obtained by replacing up with Vp in Eq. (ED . This replacement 
is needed for satisfying the charge conservation law on the disordered surface of the doped 
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III. SPIN CURRENT DUE TO MAGNETIZATION DYNAMICS 

In this section, we show spin current driven by magnetization dynamics on the disordered 
surface of the doped TI/MI junction. Here, the spin current and charge density are mutually 
related each other, because of the spin-momentum locking on the surface of the TI. 

A. Definition of spin current on the surface of topological insulators 

In order to derive the spin current on the disordered surface of the doped TI, we demon¬ 
strate the definition of the spin current. The spin current jf is defined from 

d t s a + V 8 j“ = T Q , (8) 

where s a = ^(^cr 01 ^) is the spin density, jf shows the spin current density, and T“ is the 
spin relaxation torque on the surface. Here the subscript and superscript of jf represent the 
direction of flow and spin of the spin current, respectively. From Eqs. 0-0. jf and T a 
are given by 

3 ? =y Gm(V’V) = Y^zaipe- (9) 
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with the Levi-Civita symbol e za i. In the above equation, we used the commutation relation: 

= -ihv F e za£ {V e i>^)a a + J sd ^S a cr a , 
ft/h'H] = ~ihv F e zoie a a (Vel)) - J sd S a a ar ip. 

From Eq. (J9j), the spin current is proportional to the charge density p e = where 

e < 0 is the charge of electrons. Moreover, the directions of the spin and flow of the spin 
current are perpendicular to each other because of the spin-momentum locking. The spin 
relaxation torque is derived from Eqs. AH])-©. The torque can be separated as 

T“ =7H + 73, (10) 

where T^\ and are spin relaxation torque caused by 77 ti and Usd, respectively. Here, 7^ 
and Tgd are given by 

Tri ~ ( 11 ) 

XS ( 12 ) 

We note that the definition of spin current depends on that of the spin relaxation torque 37- — . 
For example, we consider the case when the spin relaxation torque can be described by 
r" = T° + VjP“, where the polarization is an arbitrary vector with = 0, whose 

index i and a represent the direction of the polarization in the real space and that of the 
spin in the spin space, respectively. Then, the spin current can be also represented by 
J- x = jf + V (x and satisfies the conservation law as d t s a + V,;b7"“ = r a . We discuss the 
spin current defined in Eq. ©. To consider the spin current and the spin relaxation torque, 
we calculate the charge density and spin density in the following subsections. 

B. Charge density 

First, we will calculate the charge density p e in the linear response to S. p e is de¬ 
scribed by using the lesser component of the Keldysh-Green’s function, —ihG < (x,t, x,t) = 
(ip'(x,t)ip(x,t)) in the same position and time as 

p e = — ike tr [G < (a:, t, x, t)]. (13) 
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Hence, p e is given by 


p, = J]e i ( !! *-^>tr[fl 0 „(,,n)sy, (14) 

q,n 

where L 2 is the area of the disordered surface of the TI, and q = (q x , q y ) and 42 indicate 
the momentum and frequency of the localized spin S” n [y = x 1 y 1 z), respectively. Here, the 
charge-spin correlation function Ho^ is given by 

n oA<?, fi) = ^)9 k+ i a,+|] < > ( 15 ) 

k,uj 

where g k±2 w± fi is the non-perturbative Green’s function of TLt\ including 14, which is taken 
into account within the Born approximation. The retarded (advanced) Green’s function g k u 
(9k, u = [9k, J) is g iven by 

9k,u = [fux + e F - hv F & ■ (k x z) + ih/ ( 2 r)] -1 , 

where Tl/(2t) = 7miwfz/ e /2 represents the self-energy due to 14 within the Born approxima¬ 
tion. The vector A„ in Eq. (1I5|) is the vertex function, which is described by 

A 7 (q, 42) = d 7 + ^ [r(q, 42) + [f(q, 42)] 2 + • • • ] 7 £<7(. (16) 

(=0,x,y,z 

Here d 0 = 12x2 is the identity matrix and T 7 ^ is given by T 7 : 

f 7 (q, 42) = riiU i ^ ' gfc-3 a) _go , 7 gfc+3 

k 

= f 7C <r c . (17) 

The correlation function n 0l , can be decomposed into the retarded and advanced Green’s 
function by using the formula g^ u = fu,(g ki0 — g k where f u is the Fermi distribution 
function. Using the formula, we can estimate the correlation function n 0 „ on the surface of 
the doped Tl, i.e., h/(e F r) <C 1 regime as n 07 = H { r )7 + o(h/(e F r)), where is represented 
by 


n £ = £(/„+»-/„_»)? 

k,uj 

Here A™ is given by the Pauli matrix as 


A ra n a 

M"-§ 7 5i + j,„ + a- 


A~= V [i + f“ + (f“) 2 + ---] » c = J2Kc a o 


C=0 ,x,y,z 


I ™ = n,u- £ 9, 




(18) 


(19) 


( 20 ) 



Using Eqs. (THU) - (120) under the condition Qt -C 1, we can calculate in the low- 
temperature limit. Besides this, we can calculate the response function Ilg^ by postuating 
fir < 1, qi -C 1, and q 2 x = q 2 = q 2 / 2. fljj a is given by 

-O 


rrra 

1J -07 


2vr 


E -pra X ra 
1 C A iC’ 


( 21 ) 


(=0,x,y,z 

where T£ a = Y2u=o xy z T/ a ay can be expressed by 4 x 4 matrix T ra as 


( 1 — iQr — \£ 2 q 2 


pra _ 


2 ' 'iy 


2 ' 

1/92, 


V 


2 ^v 
~\hx 
0 


0 \ 

|(1 -i?LT-\t 2 q 2 ) \Pq x q y -J q x Z-£f 

\fq x q y |(1 - iQr - \t 2 q 2 ) ~{q/^ 

h i „ o h h \2 


\q X t- 

4^ J/ €pT 4 J y 6 pT v EpT ' 


( 22 ) 


epr 

, 2 _ 


°(^) 2 / 

In the above equation, we have used n^uf-KVe/(h/2r) = 1/2. Here, u e = e F /(2Trh 2 Vp) is 
the density of states at the Fermi energy on the surface of the doped TI. From the above 
equation, the magnitudes of T^ and T z ^ are negligibly smaller than that of T utl (u, p — 0, x, y ) 
for h/(e F r) <1. As a result, T y = Ysu=o,x , y \+ o(d/(e F r)) is obtained by 


1 f>2„2 ^ - 1 


rf a = 1 - iVtr - -i q d 0 + - icr a q b e abz , 


■nra 

1 [i=x,y 


1 ^1 - ittr - 


&v + —(■QatfMZ&O- 


Then, A£ a can also be estimated by using T ra as^i 

AU = [(1 - f-)- 1 ] 


7d 


(23) 

(24) 

(25) 


Therefore, from Eqs. (JT3j) and (12T1) - (12T) . the charge density p e is obtained by 


Pe = 


eVeJsdT i(ilt—q x) 


5> 

q,n 


m 


L 2 ^ q 2 ? 2 + iClr 

= -ev e J sd r£\y x ^(S 11 ^]*, 


(9y^,n - feSJn) 


(26) 


From Eq. (126|) . we find that the charge density p e is induced by d t [V x (<S^)d]z- Here 
(«S^)d is defined by the convolution of the in-plane of the localized spin and a diffusion 
propagator function T> on the disordered surface of the TI as 


{S'%{x,t) 
V(x , t) 


— J dt' J dx'V(x — x', t — t')S^(x', t'), 

_L „i(Clt—q x) j 

L 2 4^ 2 Dq 2 + iQ : 


(27) 

(28) 
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where, D = v^t/2 is a diffusion constant and (S'H)d denotes the nonlocal spin, which 
diffusively propagates by the diffusion propagator V. The V results because of nonmagnetic 
impurity scattering on the disordered surface of the doped TI. We also find that the charge 
density due to the out-of plane of the localized spin, S z , is negligible smaller than that due 
to Si 

The diffusion propagator T> satisfies the differential equation 

( d t — 2 DX7 2 )V(x — x',t — t ') = 8(x — x')5(t — t'). (29) 

We find that from Eqs. (1261) and fl27l) . the diffusive motion of the charge density obeys the 
diffusion equation: 

(d t - 2 DV 2 )p e = -ev e J sd l(y x dtS 11 )*. (30) 

The above equation means that the diffusion propagator of the charge density is caused 
by the spatial and time derivative of the localized spin, (V x d t S^) z , on the surface of the 
doped TI. When the localized spin is spatially uniform, p e is not driven by the magnetization 
dynamics. 

C. Spin current 

We will now consider the spin current due to the magnetization dynamics on the disor¬ 
dered surface of the doped TI. The spin current is proportional to the charge density [see 
Eq. (j5J)]. From the result of the charge density due to magnetization dynamics [see Eq. 
(126]) ] . the spin current is given by 

jf = ~ x ^(S 11 )^*. (31) 

This is one of the main results of this paper. From Eq. (l3Tj) . the direction of spin of the 
spin current is perfectly locked and is perpendicular to the direction of the flow of the spin 
current. The origin lies on the spin-momentum locking on the surface of the TI. The spin 
current is proportional to the coefficients, which are the density of states at Fermi energy 
v e , the s-d exchange coupling J sd , and the square of the mean-free path P. Here jf is 
proportional to the spatial and time derivative of the nonlocal spin as [V x We 

fold that the local spin does not contribute to the spin current generation. In the case 
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when the spin structure of the MI is spatially uniform, the spin current vanishes. Since the 
spin current is proportional to the charge density, we expect that the spin current can be 
arising from the accumulation of the diffusive charge density, which is given by Eq. (15UD . 
Additionally, Eq. (jUT]) indicates that the spin current is an even-function of Up, the sign of 
which depends on the helicity of electron on the surface of the TI. Therefore, the direction 
of the spin and flow of the spin current on top surface (j" top ) and that on bottom surface of 
the TI (j“ bottom ) are equal as j“ top = j“ bottom . We End that the in-plane component of the 
localized spin = S — S z z contributes to the spin current, but the out-of plane component 
of the localized spin S z z does not. We expect that its origin lies on the spin-orbit coupling 
of %t\. From (<x x p) z — a ■ (p x z) and — (z x S^) x z, the Hamiltonian 72 ti + l-i s d 
can be described by 


"Hti + 7-Lsd, = 



Vp& ■ 


(p - z x S 11 ) 

) X z 

L\ V F 

) \ 



(32) 


From the above equation, we can regard that the conduction electrons momentum p is 
shifted by the in-plane localized spin <S^: p — > p — 4^(z x «S^). The in-plane localized spin 
z x plays a role like an electromagnetic vector potential A = ^-(z x ) 13 ~ 14 . Then, the 
observable quantity should be proportional to the gauge invariant form: an effective electric 
field S = —d t A or an effective magnetic held 2? = V x A, as represented by 


S 

B 


(z x d t S H), 

evp 


—V x (z x 5 11 ). 
evp 


(33) 

(34) 


The dynamics of the in-plane component of the localized spin can be regarded as the effective 
electromagnetic held, which acts as a driving force to trigger the motion of conduction 
electrons. While, the out-of plane one S z z plays a role like magnetic helds for the conduction 
electrons and does not directly shift p in the momentum space. We expect from the difference 
of these properties of the localized spin, the contribution from S z z could be smaller than 
that from <S^. 

The spin current can be represented by using the effective electric held S and {V x 
[(£) d x z]} z = -V • (£) D as 

jf = ■ (£)n. (35) 
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From the above equation, we find that the spin current is proportional to V • (£)d stemming 
from charge density. In fact, the charge density can be represented by p e oc V ■ (£)d, as 
shown in Eq. (12611 . Here, the charge density is also proportional to V • (£)b and is similar to 
the Gauss’s law in Maxwell equations as p e = eV • E, where e is a permittivity and E is an 
applied electric field. Thus, we can interpret that the charge density and the spin current on 
the surface of the TI are generated by the divergence of the effective electric field. Equations 
m, cd. and (|35j) are the main results of this section. 

IV. CHARGE CURRENT DUE TO MAGNETIZATION DYNAMICS 

In this section, we show charge current due to magnetization dynamics on the disordered 
surface of the doped TI. Because of the spin-momentum locking, the charge current is 
proportional to the density of the spin polarization on the surface of the doped TI. We 
calculate the spin density, the charge current and the resulting spin-relaxation torque. 

A. Spin density 

To discuss the charge current, we calculate the spin density due to the magnetization 
dynamics in the linear response to the localized spin. The spin density s = ^ fri/j) is 
given by 

s “ = E e i< “-’ a ’ ) tr[rU9, n)SJ, n ], (36) 

where, n u — x,y, z ) is the spin-spin correlation function. can be calculated within 
the same formalism as in the section 3.2, and is represented by = <r M IIov From the result, 
we can obtain the spin density s. Here, s can be decomposed into two terms: s — sll + s z z, 
where = s — s z z and s z z show the in-plane and out-of plane component of the spin 
on the disordered surface of the doped TI, respectively. We find that s z z is proportional 
to d t S z , and its magnitude is negligibly smaller than that of the magnitude of gH within 
the approximation |s z |/|sll| ~ o[h/(eFr) -C 1]. Thus, the spin density can be estimated 
by s = sll + o[h/(e-pT)\ and S z does not contribute to the generation of s. The dominant 
contribution of s can also be decomposed into two terms: 

s = s L + s D , (37) 
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where s L is the local spin density due to S^, and is given by 



( 38 ) 


The local spin density s L is induced by the time-derivative of the in-plane component of the 
localized spin dt.SK On the other hand, the second term of Eq. (|37l) . s D , is the diffusive 
spin density and is given by 



(40) 


(39) 


From Eq. fl40lh s D is generated by the driving field (z x V)(V x d t (S)r>) z , which is the 
spatial gradient and the time-derivative of the nonlocal localized spin (<S^)d- The driving 
held is also described by (z x V)(V x d t (S) d) z = — V(V ■ (<S')d)]- Here, s D 


is described by the spatial gradient of the charge density, which is caused by the electron 
diffusion on the surface of the TI, as shown in Eq. (l26]h In addition, s D is also represented 


by the spin current: The charge density is proportional to the spin current, p e = A {j* — jv ), 
and s D becomes 



(41) 


From Eqs. (I38l) - (l39lh we find that s is an even-function of and is independent of the 


hclicity on the surface of the TI. Therefore, the direction of s does not depend on whether 


we are focusing on the top or bottom surface of the TI. 

We End that from Eqs. (I38l) - (l39l) . the spin is polarized not by a static magnetization 
but by magnetization dynamics. Therefore, we expect that static magnetization does not 
induce spin polarization on the surface of the TI. This seems to be anomalous property on 
the surface. The response between the spin polarization and the static magnetization on 
the surface of the TI is different from that in conventional metals: In the metals, a spin is 
polarized even by static magnetization. We will consider shortly why static magnetization 
does not generate spin polarization on the surface of the TI. The magnetization on the 
surface of the TI plays the role to shift the momentum of conduction electrons from p into 
p-^(zxSl) in momentum space. As a result, the center of Fermi sphere is also shifted 
from p = 0 into p = |^(z x S^). Then, the direction of the spin at each momentum 
are perfectly perpendicular to that of the momentum. Besides, the spin configuration in 
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momentum space does not change before and after the shift, because the direction of the 
spin at each momentum are independent of the shift. Therefore, no spin polarization is 
driven by momentum shift due to a static magnetization. 

B. Charge current 

We note that on the surface of the TI, the charge current j is proportional to the spin 
density on the surface of the TI. The charge current is represented by using the renormarized 
velocity operator as j = 2 evp(z x s)2&. From Eqs. (l38[) - fl4Q[) . the charge current is also 
decomposed into two terms: j = j h + j D as 

j L = -eu e J sd £(z x d t S 11 ), (42) 

f = eu e J sd £ 3 V[V x $<S") D ]*. (43) 

The j L is the local charge current and is induced by the time-derivative of the localized spin 

. The direction of j L is parallel to z x S H. On the other hand, j D is the diffusive charge 
current caused by impurity scattering on the disordered surface of the TI. In fact, j D can be 
represented by the spatial gradient of the charge density as j D = —2DWp e . This means that 
diffusive current is generated by the spatial gradient of the charge density on the surface of 
the TI. We note that the charge current is an odd-function of bp, so that the direction of 
the charge current on the top surface j top is opposite to that on the bottom surface jf'bottom, 
if S is same on the top and bottom surface. It is noted that, from Eqs. (1261) and (I42l) - (l43l) . 
the charge density p e and charge current j satisfy the conservation law p e + V • j 1 — 0. The 
detail is shown in Appendix [B] 

Next, we comment on the relationship between the spin current and the charge current 
on the disordered surface of the doped TI. Substituting e zai j f = vpp e /e into Eq. (1431) . we 
find that the diffusive charge current can be described by the spin current as 

= -efe 2m Vj“. (44) 

This is also the main result of this paper. We expect that the above equation displays 
the conversion between the spin current into the diffusive charge current on the disordered 
surface of the doped TI by using the spatial gradient of the spin current. The spin current 
can be converted into the diffusive charge current when the spin current depends on the 
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space on the disordered surface. The relation in Eq. (j44T) is plausible on the disordered 
surface of the doped TI, because the charge density p e is proportional to the spin current, 
and a diffusive particle current generally proportional to a spatial gradient of particles. We 
note that there is no relation between the spin current and the local charge current j L . 

C. Effective conductivity 

The charge current due to magnetization dynamics j can be also described by the effective 
electric held £: 

j = e 2 VpU e r£ + e 2 h F z/ e f*V[V ■ (£)d\- (45) 

The first term and the second term are corresponding to the local and diffusive charge cur¬ 
rent, respectively. From the above results, we will consider an effective conductivity: an 
efficiency of the charge how due to the applied effective electric held £. This is similar 
to the conventional electric conductivity: In general, the longitudinal electrical conductiv¬ 
ity is defined from dividing the charge current by an applied electric held—,. We expect 
that the current corresponds to the local current. Then, an effective longitudinal electrical 
conductivity can be defined by j L = a£, and is given by the hrst term in Eq. (1451) as 

cr = e 2 Vpi , e T. (46) 

The conductivity only depends on the parameters on the surface of the TI, and is independent 
of parameters attached to the MI. 

D. Spin relaxation torque 

We will consider the spin relaxation due to the magnetization dynamics on the surface 
of the TI. Using Eqs. (l3lf and (l37l)-(l39l). we can describe d t s a and Vjj“ as 

d t s a = ~u e J sd rd 2 S la + ^eJsJ 2 re aiz Vi[V x (47) 

= ~v e J s d 2 £odzVi[^ x 9 t (S ,|| ) D ] z . (48) 

Therefore, the spin relaxation torque T a = d t s° + V,j“ in the linear response to is 
obtained by 

T=- \ e J sd rd 2 t SW + l -u e J sd i 2 { 1 - Tdt)(z X V)[V x d t (S%] z + 0{J 2 sd ). (49) 
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The first term in Eq. (TT71) shows a local spin relaxation torque and is induced by d^S^. The 
second term is a nonlocal spin relaxation torque and is driven by (z x V)(V x dt(S^)o)z- 
The nonlocal torque vanishes when the magnetization is spatially uniform. The third term 
in Eq. (147)1 represents the higher order of J sc j and S. From Eq. (fT2]h the spin relaxation 
torque T s d is proportional to S and s, which is proportional to SK Therefore, the third 
term in Eq. (j47j) corresponds to 1~ s d within the linear response to S. We expect that the 
third term T s d — °(^ 1 ) 2 can be negligibly small and be ignored in comparison with the 
first and the second terms of T in the regime <C 1. 

The spin relaxation torque 'T is also represented by the effective electric field as 

T =^ev e v F r(d t £ x z) + ^ eu e v F £ 2 (l - rd t )(z x V)[V • (5) D ]. (50) 

The local spin relaxation is written as the time-derivative of the effective electric field and 
the diffusive one is induced by the spatial gradient of the nonlocal effective electric field. The 
electric field dependence of the spin relaxation torque on the surface of the TI is different 
from that in NM with spin-orbit interactions: The spin relaxation torque in the NM, 7 nm ; 
is proportional to the spatial gradient of the applied electric fielcPh We expect that the 
difference can be caused by the fc-dependence of the energy dispersion: The energy dispersion 
on the surface of the TI is a linear function of k , while that in the NM proportional to the 
square of k. Equations (1391) . (ITijh and (l46|) are the main results of this section. 

V. DISCUSSION 

A. Spin torque 

We will phenomenologically study the change of the magnetization dynamics in before 
and after the spin-charge generation due to the ferromagnetic resonance (FMR). Now, we 
consider a disordered surface of the MI/TI junction, as shown in Fig.l. Further, in the 
junction, a static magnetic field and ac magnetic field are additionally applied. Here ac 
magnetic field is given by a microwave irradiation, and is needed for FMR in the MI. When 
we apply this magnetic field, the magnetization dynamics is triggered in the MI, and the 
magnetization dynamics induces the spin polarization on the surface of the TI [see Eqs. 
(137)1 - (l39j) ]. Then, the induced spin polarization s plays the role of an exchange field acting 
on the magnetization in the MI. As a result, the magnetization dynamics is affected from 
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the generated spin s. Here, the mutual interaction between the magnetization dynamics 
and the induced spin are called as the feedback effect^. 

The magnetization dynamics on the disordered surface of the doped TI is described from 
the Landau-Lifshitz-Gilbert (LLG) equation as^ 

(y 

d t M = -7 n(M x H) + — M x d t M + T e , (51) 

where, M = —jf S is the magnetization of the MI, 7 is the gyromagnetic ratio, fi is a 
permeability, a is a Gilbert damping constant, H = Ho + H ac is an applied magnetic field 
on the disordered surface of the doped TI. H 0 and H ac denote a static and ac magnetic 
field, respectively. The spin torque on the disordered surface of the doped TI is given by 
Te = 2Js g a2 M x s. The torque can be decomposed into two terms: 7~ e = + 7^\ 

where 7~e = 2Js £ a M x s L and T f = 2Js g a M x s D are the local and diffusive spin torque, 
respectively. Here, a is a lattice constant on the surface of the TI. These spin torques are 
obtained from Eqs. (I37|) - (T39|) and S 11 — 35 

(52) 

Tf^( , Mx(»V)(Tx9,(Ml) D )„ (53) 

where, is the magnitude of the in-plane magnetization M H = M — M z z, and k = 
u^J^rSt/h is dimensionless coefficient proportional to J 2 sd and r. We find that oc 
M x d t M II is slightly different from the damping torque fjM x dtM ; which is a damping 
of the magnetization. We could expect that the contribution from the local spin torque 
can be observed in the experiments on the surface of FM/TI junction 27,2 -. Here, T\ plays 
the role of an anisotropic damping torque unless the static magnetic field and microwave 
are parallel to the z direction. The anisotropic damping affects the magnetic permeability: 
For example, when the static magnetic field and microwave are parallel to the ^-direction, 
then, the longitudinal magnetic permeability Xxx and y 22 are not equal each other. Here, 
' T ° seems to be a new type of spin torque 7”^ on the surface of the TI. in Eq. (153|) is 
induced by the spatial gradient of the magnetization, M x (z x V) (V x dt(M^) d) 2 . When 
the magnetization is spatially uniform, is zero and T\ is nonzero. 

Since, jf and are proportional to the charge density on the disordered surface of the 
doped TI, T° can be described by the spin current jf\ 

T? = J -^-M x [(* x V) wT]. (54) 
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The spatial gradient of the spin current induces the diffusive spin torque on the disordered 
surface of the doped TI. We expect that the contribution of jf is detected from the change 
of the half-width value, as well as the change of a shift of the magnetic resonance frequency 
through T’D [discussed in Sec JV Cj . Equations are the main results of this section. 


B. Magnetic permeability without diffusion 


Using Eqs. (HxTjh (I551h we discuss the magnetic permeability in FMR when the magneti¬ 
zation is spatially uniform on the surface of the MI/TI junction. We consider that in the 
junction, the applied static magnetic field Hq and the microwave of the ac magnetic held 
H ac are applied along the y direction: H 0 = (0, H , 0) and H ac = ( h x , 0, h z ). For an uniform 
magnetization case, the spin becomes s L ^ 0 and s D = 0, and the spin torque are ^ 0, 
but r° = 0. 

Then, the LLG equation on the surface of the MI/TI junction can be described by 

ry K II 

d t M = 7 ii(H x M) + —(M x dfM) + ^(M x $M"). (55) 

To estimate the magnetic permeability on the surface of the TI, we assume that the |iT 0 | is 
larger than the \H ac \, i.e., \ h x |<^ H and | h z |<^ H. Then, from the applied magnetic held, 
we expect that the local magnetization on the surface M = (m x , M y ,m z ) can be satished 
m x -C M y and m z -C M y . Moreover we assume that the time-dependence of the precession 
of M is given by m x oc m z oc e lQt and d t M y ~ 0. In order to solve the LLG equation, we 
take a linear approximation of 777 : rriimj r^j 0, My M{= \M\), and ~ M. Then, the 
LLG equation becomes 


d t m x = 7 n(Hm z - h z M ) + ad t m z , 
d t m z = 7 y(h x M - Hm x ) - (cc + n)d t m x , 

and the magnetic permeability is given by 


(56) 


(Xxx XzzA fh x 
'Xza: Xzz' 'h Z 

The frequency dependence of the longitudinal magnetic permeability \ xx and \ zz are de- 
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scribed by 


Xxx 

Xzz 


um(u h + iaQ.) 

(t oh + iafl) [ojh + i(a + k)Q] — XI 2 ’ 
ojh + i{a + k)H 


m .. I a r\i O 


Xxx1 


(58) 

(59) 


where, uh = 7 fiH and ujm = 7 pM are the angular frequency of the applied static magnetic 
field and that of the magnetization, respectively. The permeability Xxx and Xzz are different 
each other originating from the anisotropic spin-transfer torque 7~g. Transverse magnetic 
permeability has the relation Xxz — ~Xzx and is given by 

—iXl 


Xx 


coh + ioiXl 


Xxx 


( 60 ) 


The real part of the permeability Re[x x;c ] shows a Lorentzian profile due to the magnetic 
resonance around the resonant frequency fi r , which is proportional to H. The Im[y a . a .] 
indicates the energy absorption of the applied microwave around Xl r . Half-width value of 
Imfxzz] expresses the damping of the precessional motion of the magnetization. From Eq. 
(l58lh the half-width value AH is caused by the Gilbert damping (aM x dt.M ) and the local 
spin torque given by 


AH ~ (2a + ft)H r . (61) 

We expect that in the MI without TI, the half-width value of the magnetic permeability 
estimates AH = 2afi r . Eq. (16T1) indicates that on the surface of the doped TI, AH is 
enhanced from 2afi r into (2a + k)XI t . The origin lies on 7~g, which is triggered by the 
induced spin, where the spin is induced by the magnetization dynamics on the surface of 
the doped TI. This enhancement of AH has been verified in the recent experiment^. 

We compare AH in MI/TI junction with AH in FM/NM junction. In the FM/NM 
junction, it has been demonstrated that the enhancement of AH is triggered by the spin 
current in the NM, which is generated by the magnetization dynamics of the FM 44 -. On the 
surface of the MI/TI junction with uniform spin structure, on the other hand, spin current 
is not generated by magnetization dynamics, and the spin current does not contribute to 
AH. The enhancement of AH is caused by the spin polarization due to the magnetization 
dynamics on the surface. 
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C. Magnetic permeability with diffusion 


We will consider the surface of the TI/MI junction, where the localized spin in the MI is 
spatially inhomogeneous. The spin structure of the localized spin we consider is a spin-wave 
or longitudinal conical spin structure, which are realized in ferrimagnetic insulator yttrium 
iron garnets (YIG) or multifcrroics^, respectively. We expect that even on the surface of 
the TI, the localized spin depends on the position through proximity effects from the MI. 
Then, if we apply ac magnetic held of microwave in the junction along the y direction, we 
assume that the localized spin becomes precessional motion by the applied magnetic held. 
The localized spin on the surface of the TI can be described by 


S(x , t ) = [S cos (q ■ x — hit), S y , S sin (q ■ x — hit)], (62) 

where, S and S y are a constant coefficient independent of space. We assume | S |-c| S y \ 
and | S |~ S y . Here q = (q Xl q y ) is the momentum of the localized spin and is assumed to 
be monochromatic. The direction and the magnitude of q depends on materials of the MI. 
In the spin structure, nonlocal diffusive spin (S )d is given by Eqs. (127|) and (l62j) as 


(63) 


(S x ) D = Ag^S cos [q-x- fit]- B q ^S sin [q ■ x - Of], 

(S z ) D = Bq^S cos [q ■ x - fit] + Aq^S sin [q - x -fit]. 

The component of (S)t> is different from that of S\ (S x ) d has cos [q ■ x — fit] components, 
as well as sin [q ■ x — fit] components. Here, the coefficients A and B are obtained (see 
Appendix [All as 


4 

(Hr) 2 + (q 2 £ 2 ) 2 ’ 

(64) 

fir 

(65) 

Bq ’ n ~ (fir) 2 + (q 2 £ 2 ) 2 ' 


The coefficients A q ^ and B qt n depend on qt and fir , where parameters i and r are determined 
by the TI, and q and are chosen as characteristic values of the MI. Then, the diffusive 
spin s D is given from Eqs. (139|) as 


s D = v e J sd rf[q 2 d t (S 11 ) D - q(q ■ ^(S^d)]- 


( 66 ) 


Then, nonlocal diffusive spin torque is obtained by using (IhfH) D = —(MII/5 ^(-S , H)d as 

2 


r 


■D 


MU 


[q 2 M X ^(M'I) d - (M X q)(q ■ d t (M") 


D 


(67) 
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The first term of Eq. (16?|) plays the role of an anisotropic damping torque, which is similar 
to 7^ l . The direction of this torque is independent of the direction of q. On the other hand, 
the second term of Eq. (jU71) is proportional to (M x q)(q ■ d t (M^)- D ), and its direction 
depends on q. 

Using Eqs. (1S21)- (R)T1) . we consider the magnetic permeability affected by T’^.Then, the 
LLG equation is given within the linear order of the magnetization as 


d t m x = uj H m z - h z uj M + oid t m z 

d t m z = h x uj M - u H m x - (a + n)d t m x + nq 2 £ 2 d t (M^ x ) D . 


The last term of the above equation is caused by 7^°. In order to discuss the permeability due 
to 7)?, we consider when the momentum has the y-component (q = qy), whose direction 
is parallel to the applied static magnetic field. That means that the spin structure we 
consider is a longitudinal conical spin structure. From the above equation, using = 

Af h Qrn x — B q ^}rri z , we obtain the magnetic permeability as 



Here the longitudinal and transverse magnetic permeability are given by 


X xx (<lM) 




Xxz(<l > 


{ojh + i<yQ)coM 

[uj h + «afi] [oj h + i(a + /c g) n)fi] - C?,o ^ 2 ’ 
c o H + + k q p)Q D 

u H + iaXl Xxx ’ 

CqftXzx' 


( 68 ) 


(69) 

(70) 

(71) 


The obtained permeability is different from that in Eqs. (l58lb()60lh The difference is caused 
by coefficients k (h Q and 

K,.n = K(l (72) 

(q,n = 1 + nq 2 £ 2 B qM . (73) 


The K q) n and Q Ch n depend on qt and Hr. If q — 0, one can demonstrate k q ,a = k and 
= 1) and the magnetic permeability in Eqs. (j69l) - (17T]) equal to that in Eqs. (I58 l) -(l60 lh 
respectively. Figure 2(a) indicates the Hr dependence oik/n for several momentum q£. The 
parameter k/k approaches to 0 from k/k = 1 with increasing Hr: In the case for k/k — 0, 
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qi 

FIG. 2: (Color online) (a) The fir dependence of k/k for several qi. (b) The q£ dependence of the 
£ for several frequency (2nf = 1, 2, 5,10 GHz) in the fixed relaxation time (r = 0.087 ps). 

Tp and TP are canceled out each other, and the spin torque T e vanishes. On the other 
hand, k/k = 1 means that 7p is zero and 7p is nozero. The relation k/k significantly 
changed when ( q£ ) 2 ~ Or is satished. Figure 2(b) shows the q£ dependence of ( for several 
angular frequency of the applied ac magnetic held O, where we take a realistic relaxation 
time (r = 0.087 ps —). The magnitude of ( changes around qi ~ vTFF and approaches to 
( = 1 with increasing qi. 

Figure 3 (a) and (b) show the Or dependence of the real and imaginary part of the 
longitudinal magnetic permeability, Re[yJ?J and — Im[xJ?J, respectively. In Figs. 3(a) and 
(b), we choose realistic parameters of the TI: l = 40nm, vp = 4.6 x 10 5 m/s^, r = 0.087 ps, 
£ = 0.003, and v-p/v-p = 0.997. Besides, we choose material parameters of the ferromagnets^: 
a = 0.015, J S( i ~ 6meV, ~ 0.3, and Fermi wavenumber kp = 3.9 x lO 8 !!!” 1 . Then, 
n = 0.008 is obtained. The frequency of the magnetization u>m/{ 2tt) — Im — 0.28 GHz 
in the MI. The magnitude of the frequency is evaluated by the material parameters of the 
permalloy 28 . 
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FIG. 3: (Color online) (a)-(b) The frequency dependence of the real part and imaginary part of the 
longitudinal permeability, Re[xJ? x ] and — Im[y xa ,], for a fixed Gilbert damping constant (a = 0.015), 
the relaxation time [r = 0.087ps], the anisotropic damping constant [k = 0.008], and the frequency 
[fn = 1.6 GHz and /m = 0.28 GHz] for several momentum q£. 

We plot the Re[y xx ] and —Im[y; xx ] functions as the frequency of the ac magnetic field for 
several momentum of the localized spin, when we take the frequency of the static magnetic 
held uj h /(2tt) — f H — 1.6 GHz. The magnitude of fir can be estimated as fir ~ uj h t = 
8.6 x ICC 3 . For q£ < 0.03, the permeability Re[y x J and — Im[y x J are not dramatically 
changed from that without diffusion, Re[y xx ] and — Im[xxx], respectively. The reason is due 
to the prohle of k/k and (: Both k/k and ( are about 1 within q£ < 0.03 in f It ~ 8.6 x 10” 3 . 
While qt is near q£ = 0.03, Xxx changes: The magnitude of Xxx increases from that of Xxx- 
Besides, the resonant frequency of Xxx increases from that of Xxx [ see Q? — 0.03 in Figs. 3 
(a)-(b)]. The change of the resonant frequency is also shown in Fig. 4 (a). We find that 
when (qt) 2 ~ fir is satisfied, k/n and ( deviate from 1, as shown in Figs. 2(a)-(b). After 
increasing ql from qd. = 0.03, the shifted resonant frequency gets back again. The magnitude 
of the permeability increases with increasing q£. 

The half-width value is also changed by the magnitude of Inr[x^ x ] for several momentum 
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0 0.02 0.04 0.06 0.08 0.1 

FIG. 4: (Color online) The qt dependence of the resonant frequency (a) and the half-width value (b) 
normalized by the angular frequency of the applied magnetic field for several frequencies (27 t/h =1, 
2, 5, and 10 GHz) in 2 a + k = 0.038 and 2 a = 0.030. 

qt. This trade off between the half-width value and qt for several frequency (27 rfn = 1, 
2, 5, and 10 GHz ) is shown in Fig. 4(b). In the case for r = 0.087 ps, the half-width 
value Afi significantly changed around the qt , which satisfies qt VQt (e.g., qt = 0.03 in 
u>h = 10 GHz). Figure 4(a) indicates qt dependence of the resonant frequency rate Q r /uH- 
The VL r /ujH changed from Q. t /ujh ~ 1 into fl r /cuH = 0.998 around the qt, which satisfies 
qt vThr. 

We discuss the role of diffusive spin torque from Figs. 4(a) and (b). Figure 4(a) 
shows that the resonant angular frequency Q r tends to decrease with increasing qt. The 
decrease of can be caused by the increase of and the change of ( q ^ is caused by T”j?. 
Therefore, we expect that 7"^? plays a role as a field-like torque to shift the resonant fre¬ 
quency. Then, Fig.4 (b) indicates that the half-width value tends to decrease with increasing 
in Eq. (73), which is caused by 7~J?. It means that the damping of the magnetization 
dynamics is reduced by on the disordered surface of the doped TI. Thus, behaves 
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TABLE I: A brief summary of the role of the spin torques and the half-width value in several 
spin structures on the surface of the TI, when the direction of Hq and the propagation of H ac are 
parallel to y-axis. 


Uniform 

Transverse conical ( q y 

= 0) Longitudinal conical ( q y / 0) 

7^ l Damping torque 

Damping torque 

Damping torque 

T D 

- 

Damping torque and field-like torque 

AQ/Q r 2 a + k 

2 a + k 

2 a + h 


both the damping torque and the field-like torque. The role of the spin torque in several 
spin structures are summarized in the Table I. From the table I, we can expect to distinguish 
the half-width value contributed from 7^ L and by tuning of the magnitude of the applied 
magnetic held H 0 . The reason is why the longitudinal spin structure changes an spatial 
uniform ferromagnetic structure if when we apply a strong magnetic held, which broke the 
longitudinal spin structure. Then, qi of the spatial uniform spin structure can be regarded 
as ql = 0. Therefore, we expect that if when we apply the strong applied magnetic held 
in the longitudinal spin structure, the half-width value with a finite qi changes into the 
half-width value with qi = 0 as AL2/L2 r . = 2a + k —» 2a + k [see Fig.4 (b)]. 

For magnetization dynamics due to the magnetic resonance, we need to apply magnetic 
holds in the MI/T1 junction. Then, the spin-charge generation and transport are triggered 
not only by the magnetization dynamics, but also by the applied magnetic held. We will 
estimate when the contribution due to the magnetic held can be relevant. The contribution 
from the applied magnetic held can be described by the Zeeman effect, "Hz = — 2hq f dxBs , 
where B is the applied magnetic held and couples with conduction electrons spin on the 
surface of the TI. The contribution from B can be treated within the same formalism in 
sections 3 and 4 by replacing S —* S + 2(hr) / J s< j)B in Eq. (J3J). As a result, the spin-charge 
generation and transport due to "Hz and 1-L s d are obtained by replacing S — > S + 2 (fry/ J s d)B 
in Eqs. (j26|) . (jHTlh (I38l) -(l39l). and (j42l) - (l43lh We expect that the contribution from "Hz can be 
ignored compared with that from T-Lsd, when the energy scale of the Zeeman effect is smaller 
than that of the exchange energy on the surface of the TI as 2h^\B\/(J sd \S\) 1. The 

\B\/J sd value can be estimated by \B\/J sd <C l/(2h/y) ~ 4x 10 4 T-eV _1 in S ~ 1. Then, the 
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magnitude of the exchange coupling J s d can be estimated. Using realistic parameters, the 
mean-free path £ = 40nm, Fermi velocity vp = 4.6 x 10 5 m/s^, we obtain the h/r ~ 15meV 
and J s d 15 meV, which is requested for the perturbation condition J sd r/h 1. Then, if 
\B\ 10T is satisfied, we need to consider the contribution from "Hz- 

At the end of this section, we estimate the spin density due to the dynamics of the 
longitudinal spin structure S = (sin 6 cos fit, 1, sin 9 sin fit) in the case of 6 <C 1 around 
resonance angular frequency fi ~ 1 x 10 10 s -1 . The magnitude of the spin depends on the 
regime of q 2 £ 2 fir or q 2 £ 2 fir. In q 2 £ 2 fir regime, the magnitude of the nonlocal spin 
can be negligible small compared with that of the local term, and the spin is estimated by 
|s| ~ 1.6x 10 -10 A -2 at 9 rs./ 0.1 rad in the FMR. Then, we find that the magnitude of the spin 
due to the spin-pumping is smaller than that of the spin due to the applied electric fields. 
On the other hand, in q 2 £ 2 fir regime, the local and nonlocal spin vanishes each other 
even in the presence of FMR. From the results, we expect that the change of the magnitude 
of the spin dependent on q 2 £ 2 /fir can be measurable for several applied magnetic fields, 
because the inhomogenous spin structure (q£ ^ 0) changes into an uniform spin structure 
(q£ r\_/ 0) by using an applied strong magnetic held. 

D. Spin current and charge current 

We will discuss the spin current on the surface of the disordered MI/TI junction compared 
with that in the FM/NM junction. The spin current due to the spin-pumping j“ FM / NM in 
the FM/NM junction is triggered by the magnetization dynamics a&h- 1 ^ 

= b\7AS a + 0(S 2 ), (74) 

where b is a coefficient dependent on materials. It is similar to the spin current in Eq. (1741) , 
that the spin current is proportional to a time-dependent magnetization and the spin current 
vanishes when the magnetization is spatially uniform. The direction of the spin (a) and the 
how ( i ) of the spin current in Eq. (1741) are not related each other. On the other hand, spin 
current on the surface of the TI, whose direction of spin and how are perfectly perpendicular 
to each other. The difference lies on the spin-orbit interaction. The jf in Eq. (13TT) includes 
the contribution of the spin-orbit interaction, which is absent in Eq. (T74ji does not. 

Charge current due to the spin-pumping in the FM/NM junction is also given by the 
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magnetization dynamics and Rashba type spin-orbit interactions. When the magnetization 
is spatially uniform, the charge current j"fm/nm becomes^ 

Jfm/nm = a x (S x d t S ), (75) 

where a is a constant vector including the contribution from the spin-orbit interaction. The 
charge current in Eq. (175|) is proportional to the second-order of the localized spin S. That is 
different from the charge current on the surface of TIs. The charge current on the surface of 
the TI is proportional to the localized spin d t S", as shown in Eq. (I43H . The difference is due 
to the property of the localized spin: the localized spin plays the role of the effective vector 
potential on the surface of the TI. We note that the frequency dependence of these charge 
current is also different; j"ti oc dtS^ oscillates with time of the localized spin in the FMR, 
but j’fm/nm oc Sxd t S does not. For example, the ac current j'ti oc (cos fit, sin fit, 0) is given 
when we apply the magnetic held parallel to the z direction on the TI/MI junction. On 
the other hand, the dc current Jfm/nm oc fl(0, 0,1) is obtained when we apply the magnetic 
held parallel to the z direction in the NM/FM junction. 


VI. SUMMARY 

We have studied the spin-charge generation and transport due to the magnetization dy¬ 
namics on the disordered surface of the doped TI/MI junction. The spin current jf i is 
proportional to the charge density p e and the direction of its spin and its how are perfectly 
perpendicular to each other, because of the spin-momentum locking on the surface of the 
TI. We have found that jf and p e are induced by the time- and spatial-dependent of nonlo¬ 
cal magnetization dynamics, which is affected by nonmagnetic impurity scatterings on the 
disordered surface of the doped TI. These results of jf and p e are shown in Eqs. (l3Tj) and 
f[26jl . respectively, jf and p e is induced except when the magnetization dynamics is spatially 
uniform. We have also shown the induced spin s and charge current density j due to the 
magnetization dynamics. Because of the spin-momentum locking, the spin s and charge 
current density j are proportional to each other. The s and j are generated not only by 
the local magnetization dynamics, but also by the nonlocal magnetization dynamics with 
the diffusion on the disordered surface of the doped TI. These results of s and j are shown 
in Eqs. (I38jl - (l39l) and (I42l) - (l43l) . respectively. A brief summary of the local and nonlocal p e , 
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TABLE II: A brief summary of the charge, charge current, spin, and spin current density due to 
the magnetization dynamics on the disordered surface of the TI. These terms are driven by the 
effective electric field S. 


Local term 
Nonlocal term 
Driving field 


charge density p e current density ji spin density s a 
z x dtS^ dtS 11 

[V x d t (S^) D ] z V[V x dt(SW)T>] z (zxV)[Vx9 i (S») D ] z 
V • (£) d e, V 2 (5)d zxs, (zx V)[V-(5)d] 


spin current density jf 


[V x 5 t (SH) D ] 2 

V • (5) d 


jf, s, and j due to the magnetization dynamics are represented in Table. II. As a result, 
we have discussed the modification of the magnetization dynamics before and after these 
spin-charge generation and transport on the disordered surface of the doped TI. These spin- 
charge generation and transport can be detected from the half-width value of the magnetic 
permeability in the magnetic resonance in the MI/TI junction, as discussed in section V. 
The magnitude of a Gilbert damping constant a in ferromagnetic insulator is smaller than 
that in ferromagnetic metals. Then, we can easily detect the change of the fn dependence 
of the resonant frequency and half-width value, which are shown in Figs. 4 (a)-(b). 

The preparation of the hybrid system with the ferromagnetic insulator deposited on the 
surface of the TI, EuS/B^Ses, has been reported^, where the magnetic moment of Eu 
locates at the interface between the EuS and B^Ses. If the magnetization dynamics of 
the magnetic moment of Eu is triggered by an applied magnetic held, the spin density 
and charge current can be induced on the surface of the TI. Additionally, the magnetic 
distribution of Eu has a magnetic domain, which is spatially dependent on the position on 
the surface. Therefore, when we move the magnetic domain by using an applied magnetic 
held, the charge density and the spin current are also triggered only around the magnetic 
domain. Recently, magnetic insulator with noncoplanar spin structure has been reported. 
For example, magnetoelectric insulator Cu 20 Se 03 has spatially dependent spin structure, 
and is called skyrmion, which is topologically protected magnetic spin vortex-like object^ 1 ^. 
If one can prepare the vortex-like spin structure deposited on the surface of the TI and can 
trigger the magnetization dynamics of the skyrmion, we expect that the charge and spin 







currents are driven by the magnetization dynamics of the skyrmion. Moreover, the spatial 
distributions of the charge density and the magnitude of the spin current depends on the 
position of the skyrmion, because the spatial derivative of the localized spin depends on the 
positions in the skyrmion. We comment that the induced spin and charge currents could be 
independent of the polarity of the skyrmion, since these current are triggered by the in-plane 
component of the localized spin [see Eqs. (T3Tj) and (l42]) - (l43j) ]. Then, we expect that in the 
Mis deposited on the surface of the TI, the magnetization dynamics induces not only the 
local spin-charge generation and transport, but also the diffusive one. Our obtained results 
will enable the applications of TI nanomembrane in spintronics devices. 
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Appendix A: Derivation of Eqs. (1651) - (1661) 

We show details of calculation of coefficient A q ^ and B q ^ of nonlocal localized spin (S')d- 
The (S) D is given by 



To substitute S = S( 1, 0, —7)e* ( ^ <3, £C oc e^ 9 ' 33 nt ) into the above equation, we calculate 



gi(q-x—Lot) 


(A2) 


q 2 £ 2 — iflr 
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The resulting ( S x )n is obtained from the real part of S(n )d in the above equation. In the 
same way, (S z )d is obtained from the real part of —iS(n) d- Thus, from the above equation, 
the coefficient of A q Q and B q ci are derived as Eqs. ( 1551 ) - ( 1551 ) . 


Appendix B: Charge conservation 


To check the validity of the spin current and the charge current we calculate, we use the 
charge conservation law d t p e + V • j = 0. The charge density d t p e is given by 


d t pe = 


er'g J sy/t 

L 1 2 3 




[Qt-q-x] 


i£Q 2 


L qfo 


q 2 £ 2 + iCtr 


- g x S’_ n ) 
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where £ = vpr is the mean-free path. The resulting V ■ j becomes 
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= -a,/,,. (B2) 


Therefore, the p e and j satisfy d t p e + V • j — 0. 
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